Nonlinear coherent loss for generating non-classical states 
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Here we discuss generation of non-classical states of bosonic mode with the help of artificially 
designed loss, namely the nonlinear coherent loss. We show how to generate superpositions of Fock 
states, and how it is possible to "comb" the initial states leaving only states with certain properties 
in the resulting superposition (for example, a generation of a superposition of Fock states with 
odd number of particles). We discuss purity of generated states and estimate maximal achievable 
generation fidelity. 
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I. INTRODUCTION 

Artificially designed nonlinear loss is a rather novel 
and highly promising way to generate deterministically 
non-classical states. Nonlinear loss possesses a number 
of advantages over more usual nonlinear unitary dynam- 
ics. First of all, it can turn a mixed initial state into a 
final pure one, which is, generally, impossible with uni- 
tary transformations. Then, the artificial loss is, actually, 
turning a greatest enemy of the non-classical state gener- 
ators into their biggest ally. In schemes of generation of 
non-classical states via unitary dynamics linear losses are 
prone to destroy quantum features rather quickly (here 
one might recall famous " Schroedinger-cat states" turn- 
ing into mixtures with the rate proportional to the aver- 
age particle number of the state However, nonlinear 
losses can force the system into a certain non-vacuum 
state ("pointer state" 0,0]), and make it remain there 
despite possible perturbations (which was epitomized in 
a concept of "protecting" quantum states [3t9)- Then, 
the nonlinear loss generation scheme can be extremely 
robust, and is able to turn a wide variety of initial states 
into a predefined output state Q. For example, one can 
produce " Schroedinger-cat states" [1] , and Fock states Q 
this way. The nonlinear loss can be a flexible and efficient 
tool of quantum computation p^ . 

Recently, significant progress was achieved in develop- 
ment of practical realizable schemes of artificial nonlinear 
loss. It was shown that it is possible to produce a wide 
variety of nonlinear losses for vibrational states of ions in 
magnetic traps [ill, [13 , or atoms in optical lattices [11] . 
Using a combination of optical nonlinearities, one can 
produce nonlinear losses of electromagnetic field modes 
[l4| - [l7| . Implementing correlated linear loss, one can also 
produce an effective nonlinear loss in Bose-Einstein con- 
densates [3 ■ Correlated loss gives an opportunity to de- 
sign nonlinear losses for the generation of single-photon 
states and st rong ly sub-Poissonian states in multi-core 
optical fibers [l9| . 

This work is devoted to analysis of possible states 
which can be produced with the help of designed non- 
linear coherent loss quite commonly encountered in the 
schemes mentioned above. This kind of loss is able to 
produce a wide range of non-classical states. We discuss 



conditions and prerequisites for generation of such states, 
provide an estimation of fidelity and purity of generated 
states. In particular, we show how to generate an ar- 
bitrary Fock state with fidelity arbitrarily close to unity 
using nonlinear coherent loss and an input coherent state. 
Also, we analyze the generation of finite or countable su- 
perpositions of Fock states. It appears that one cannot 
reach a unit fidelity for generation of superposition of two 
Fock states. The fidelity of such superposition generation 
gets worse when the states of the pair are strongly differ- 
ent in particle numbers. On the other hand, we show that 
it is possible to cut a countable set of Fock-state compo- 
nents from the initial state with almost perfect fidelity, 
thus effectively " combing" the input state and producing 
almost pure state as the result. 

The outline of the paper is as follows. In Section [II] 
we describe nonlinear coherent states and introduce a 
nonlinear coherent loss (NCL). Here also the the master 
equation for coherent losses is given. 



In Section iHl] we investigate general properties of the 
NCL and analyze their connection with features of a sta- 
tionary state which is to be obtained with the help of the 
considered kind of the NCL. 



In Section IIVI we discuss the generation of pure Fock 
states and consider the problem of coherence preserving 
for two important classes of states that can be generated 
by the discussed approach starting from a coherent state: 
superpositions of two Fock states (either including vac- 
uum state or not) and superpositions of states with pho- 
ton numbers distributed with equal intervals (e.g. states 
with even or odd photon numbers). 

In Section |V] the conditions, necessary for complete 
preservation of coherence, are derived. It is shown that 
coherent states taken as starting states for the considered 
type of evolution do not satisfy the conditions exactly 
and, therefore, cannot be used for creation of completely 
pure nonclassical states by nonlinear absorption. How- 
ever, an optimal classical starting state belongs to the 
class of coherent states, and the final state, arbitrarily 
close to pure state, can be generated by choosing suffi- 
ciently high amplitude of the starting coherent state. 
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II. NONLINEAR COHERENT LOSS 

In our work we consider dynamics of a single bosonic 
mode, described by the following standard master equa- 
tion in the Lindblad form: 



dp{t) 
dt 



7 (2Lp{t)L+ - L+Lp{t) - p{t)L+L 



(1) 



where 7 > is the decay rate and L is the Lindblad 
operator. We assume that our bosonic mode is described 
by the creation and annihilation operators and a. We 
shall consider the non-unitary dynamics generated by the 
following general class of Lindblad operators 



L = af{h) 



(2) 



where n = a+a, and f{n) is a non-negative function. 

The Lindblad operator 1^ can be considered as the 
annihilation operator of so-called /-deformed harmonic 
oscillator with the commutation relations 

[L,L^]^{h+l){f{n + l)f-n{f{n))\ 

Eigenstates of the operator L were termed "nonlinear 
coherent states" [1^ (it is curios to note that a specific 
subclass termed " Mittag-Lefler coherent states" does ac- 
tually arise in micromasers in presence of loss and inco- 
herent pump [24] ). For that reason, we refer to the de- 
cay described by the Lindblad operator ([5]) as " nonlinear 
coherent loss" (NCL) in further consideration. It is inter- 
esting that any pure state non-orthogonal to an arbitrary 
Fock states can be exactly represented as a nonlinear co- 
herent state [H. If it is orthogonal to some Fock states 
(for example, if this pure state is a finite superposition of 
Fock states), then one can still build a nonlinear coherent 
state closely approximating the state in question l51]. 

Thus, NCL looks highly promising for non-classical 
state generation. In Ref. 5] is was shown how to build a 
function / leading to an approximate generation of an ar- 
bitrary Fock state from the initial coherent one. In more 
recent work it was shown that the single-particle 
Fock state, can be generated with arbitrarily high 
fidelity from the initial coherent state by the NCL with 
the Lindblad operator 1^ with f{n) = h — 1. 

It should be emphasized that NCL design is completely 
realistic and could be realized in practice. For vibra- 
tional states of ions in magnetic traps NCL is already 
realized [Til Ull • Recently it was shown how to produce 
the NCL by simple adjustment of well parameters in the 
three-well trap configuration for Bose-Einstein conden- 
sates [18] . Very recently the new realistic way of realiz- 
ing NCL in multi-core optical fiber was suggested 
It was shown that for experimentally realistic values of 
Kerr nonlinearity of chalcogenide fibers with subwave- 
length core (which is 10^ times higher than Kerr nonlin- 
earity of a conventional fused silica optical fiber [20, 21]), 
it is possible to realize NCL in such a scheme and achieve 
a deterministic generation of a single-photon state [l9| . 



Note that the losses and imperfections of the scheme, 
generally, do not spoil the desired form of NCL (however, 
they can lead to the appearance of other losses, both 
linear and nonlinear, spoiling the effect of NCL). The 
form of NCL is defined by the nonlinearity present in the 
scheme. For illustration in the Appendix A an example 
of the nonlinear loss appearance (and NCL appearance, 
in particular) is given for the system of nonlinear bosonic 
modes coupled to the strongly dissipative mode. 

In this work we do not intend to discuss practical real- 
ization of NCL in more details. Our aim is more funda- 
mental; we want to discuss the very possibilities offered 
by NCL. Further in this work we discuss general limita- 
tions on the states that can be generated by NCL starting 
from a coherent state. We demonstrate when it is pos- 
sible to generate pure superpositions of Fock states, and 
when losses lead to the decrease of the states purity. 



III. GENERAL PROPERTIES OF FUNCTION 
/(n) AND CORRESPONDING STATIONARY 
STATES 

In order to analyze dependence between general prop- 
erties of the function f{n), describing nonlinear absorp- 
tion, and stationary states that can be obtained as a 
result of evolution, characterized by the master equation 
Eq. ([l} , it is convenient to decompose the density matrix 
in terms of Fock states: 



(3) 



Then, the master equation Eq. ((T]) leads to the following 
system of equations for density matrix elements: 



dt 



= 2jF{n + l)F{m + l)p„+i,™+i(t) - 
-j{F\n)+F\m)}pn,n{t), 



(4) 



where F{n) = y/nf{n) > 0, F(0) = 0. 

Note that, according to the system (g]), different diag- 
onals of the density matrix evolve independently. This 
fact can be made more apparent by introducing notation 



Pn,n+k{t) = Ckik{n,t), 



(5) 



where constants Ck can be arbitrary and will be fixed 
later. Quantity ^fc(n,t) satisfies the following equation: 



dt 



2-fF{n + l)F(n + k+ l)^k{n + l,t)~ 
-7{^'(n)+F2(n + fc)Kfc(n,i), 



(6) 



which does not contain £^k'{n',t) for k' 7^ k. 

The density matrix is the Hermitian one. It is com- 
pletely defined by elements with m > n. Further, 
we will take into consideration only the main diagonal 
of the density matrix and the diagonals, lying below the 
main diagonal. We assume fc > in Eqs. (O, ([6]) and 
similar equations. 
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Lemma 1. // quantities ^k{n,t) are positive (non- 
negative) for all n at t = Q, they remain positive (non- 
negative) for all t > 0. 



Proof. Eq. ([B]) can be rewritten as 

dik{n,t) 



dt 



2-/F{n + l)F{n + k + + l,t)x 

^jt{F^{n)+F^{n+k)-F^{n+l)-F'^{n+k+l)) 



(7) 



where |fc(n,t) = ^^(n, i)eT*(^'(")+^'("+'=». Right-hand 
side of Eq. ([7]) is non-negative for positive or non- 
negative vahies of quantities ^k{n,t). It leads to non- 
negativity of derivatives on the left-hand side of Eq. ([7]) 
and non-decreasing character of evolution of the quanti- 
ties if they are initially non-negative. The expo- 
nential factor, connecting ^fc(n,t) and ^k{n,t) is strictly 
positive, and, therefore, quantities £,k{n,t) (together with 
preserve their positivity (non-negativity) during 
evolution. □ 

If the initial state of the considered field mode is the 
coherent state |a), the condition of Lemma [T] can be sat- 
isfied by setting Cfc = (a*)'' in Eq. Then 



Un,0) = Pn.n+k{0)/ck = |ap"e-l"l /Vn!(n + fc)! > 
for all n, k. 

It should be noted that any density matrix can be rep- 
resented in diagonal form in terms of coherent states us- 
ing the Glauber-Sudarshan P-representation ^ 



p{0)= / P{a)\a){a\, 



(8) 



where P{x) is the Glauber-Sudarshan quasiprobability 
distribution. Moreover, one can closely approximate the 
state in question using a discrete set of coherent-state 
projectors on a square lattice [l^. Therefore, one may 
conclude: all results derived here for an initial coherent 
state conditions will be valid for an arbitrary initial state. 

Lemma 2. //^^(n, 0) > for all n, for the stationary 
value ^k{ni) = limt_j.oo ^kini,t) of the quantity S,k{ni,t) 
to be non-zero, it is necessary that F{ni) — F{ni + k). 



Proof. Summing Eq. 
relation: 

d ^ 



over n leads to the following 



F{n-^k)yUn,t), 



(9) 

where we have taken into account that F{0) = 0. For 
the stationary state the left-hand side of Eq. ^ equals 
zero. According to Lemma [U each term on the right- 
hand side of Eq. ^ is non-negative, and the sum can be 
equal to zero in stationary state, when for each n either 
$fc(n) = limt^oo ^(n, t) = or F{n) = Fin + k). □ 



Theorem 1. Density matrix element Pn,n+k{t) can have 
non-zero value limt_j.oo Pn,n+fe(i) ^ m the station- 
ary state of evolution, described by the master equation 
Eq. (gp, only if 



F{n) = F{n + k) = Q. 



(10) 



Proof. If the initial state of the field mode is a coher- 
ent state p(0) = |a)(a|, conditions of Lemmas [TJ [5] 
are satisfied by setting Ck = {a*)^. Therefore, for 
limi_+oo Pn,n+k{t) 7^ it is necessary that F{n) = F{n + 
k). Eq. ([6]) implies that the following relation is satisfied: 



n-l 



dt 



J2 (m, t) = 2-/F{n)F{n + k)(k (n, t) 







n-l 

7 E {^M - + k)f ik{m, t). 

m=0 



(11) 



In the limit t — > oo Eq. (fTTI) is transformed into 

2^F^(n)Un) = Q. (12) 

because, according to Lemma [21 for all m either F{m) — 
F{m -I- fc) = or 5fe(rn) = limt^oo £,k{m,t) = 0. There- 
fore, for limt^oo Pn,n+k{t) = Ck£,k{ri) ^ 0, it is necessary 
that F{n) = F{n + k) and F{n) = 0. 

Taking into account that any initial state of the con- 
sidered mode can be represented using P-representation, 
one concludes that the quantity lim4_j.oo Pn,n+fe(i) can 
have non-zero value, only if it is non-zero in the sta- 
tionary state for at least one coherent state \a) taken as 
the initial state. Then the first paragraph of the Proof 
is applicable, and we obtain F{n) — F{n + k) = as 
necessary condition. □ 

In order to derive explicit expressions for stationary 
values of non-zero elements of the density matrix, one 
can rearrange system of equations ([6]) in the following 
way. By introducing quantities 



Tk{n) 



2F{n)F{n -f- k) 



F^{n) + F'^{n + k) 
Eq. ([B]) can be transformed as 
d£,k{n,t) 



dt 



^Tk{n + l)$fc(n + l)6(ri + l,t)- 
-7*fc(rt)a(n,i), 



(13) 



(14) 



where $fe(n) = F"^ (n) -\- F'^ {n -\- k) . Then, it is quite easy 
to show that the following equality holds: 

d 

^ E ^k{m,t)Tk{ni + 1) • ... • Tk{m - l)Tk{m) = 

m— Til 

= 2-fF{n2)F{n2 + k)^k{n2,t) - 7$fe(ni)efc(ni, t). 

(15) 

Let the number ni correspond to the quantity 
^fc (jii , t) , which has non-zero stationary value (ni ) = 
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limt^ao £,k{ni,t) (i.e. the conditions F{ni) = and 
F(ni + fc) = are satisfied). We can choose 77,2 to be 
the minimal number greater than ni, for which at least 
one of the following equations is satisfied: F(n2) — 
or F{n2 + k) — (712 can be equal to infinity; then 
£,k{n2,t) — > 0). For this choice of ni and n2 the right- 
hand side of Eq. dTSl) equals zero. The left-hand side of 
Eq. (jl5p remains constant during evolution. In the limit 
i — > 00 we obtain 

^kM = a(ni, 0) + a(ni + 1, 0)Tfe(ni + 1) + 

+ ...+a(«2-l,0)rfe(ni + l)-... •rfe(n2-l). ^ ' 

Using Eq. ([5]) and returning to the density matrix ele- 
ments, one can rewrite Eq. (jl6p as 

Pnuni+k{<X)) = Pnuni+k{0) + 

(0)Tfc(ni + l) + ...+ (17) 
+ Pn2-i.n2+k-iiO)Tk{ni + 1) • ... • Tk{n2 - 1). 

The obtained expression for non-zero elements of the 
stationary density matrix can be interpreted in the fol- 
lowing simple way. The master equation ([1} in the form 
(HI describes "flow" of amplitudes of density matrix ele- 
ments along diagonals in the direction of photon number 
decreasing. The " transmittance" of the transition be- 
tween pn,n+k and pn-i,n+k-i equals Tk(n) (see Eq. (0). 
Elements Pn,n+k with F{n) = F{n+k) = "accumulate" 
the flow (i.e. they do not transmit it to the next elements 
Pn-i,n+k-i)- Elements Pn,n+k with either F{n) = or 
F{n + fc) = (but without the two conditions being 
satisfied simultaneously) neither transmit the flow, nor 
accumulate it; for such elements Tk{n) — 0. This inter- 
pretation is illustrated in Fig. [1] 

It should be noted, that, according to the definition 
Eq. ([T^ of " transmittances" Tk{n), amplitudes of diago- 
nal elements of density matrix are transmitted perfectly: 
To{n) = 1 (i.e. the trace is expectedly preserved). In 
order to attain maximal coherence, transmittances for 
non-diagonal elements must be also equal to unity (at 
least those present in expressions of the form of Eq. PT)) 
for non-zero elements of the density matrix). Therefore, 
the condition of preserving maximal coherence in the sta- 
tionary state is 

F{n)^F{n + k), (18) 

for density matrix elements giving non-zero contributions 
to non-zero elements of the stationary density matrix. 



IV. EXAMPLES 

A. Generation of Fock states 

The most "natural" nonclassical states that can be 
generated by NCL with arbitrarily high fidelity are pure 
Fock states. Here and further in this Section we assume 



(a) 

Tk{ni + l) Tk{ni + 2) ri.(77 2 + 1) = 




FIG. 1: (a) Interpretation of system evolution in terms of 
"flow" of amplitude of density matrix elements. In the shown 
example F{n\) = F{ni + k) = and "flow" is accumu- 
lated by the density matrix element pn^^.m+k', f (1-2) = 
and F{n2 -I- fc) / or ^(712) / and F{n2 + k) = and 
"flow" is "absorbed" at the transition between pn2,n2+k and 
Pn2-i,n2+fc~i without accumulatlou. (b) Schematic repre- 
sentation of the density matrix for the case, when F{0) = 
F{ni) — and F{n) 7^ for n 7^ 0, ni. Black squares repre- 
sent density matrix elements, which are non-zero in the sta- 
tionary state. Grey squares represent elements, which give 
contribute to non-zero elements of the stationary density ma- 
trix. Arrows indicate direction of "flow". Thick lines corre- 
spond to transitions with zero transmittances. 



that the initial state is the coherent state, \a), with the 
amplitude a. 

Suppose that the function f{n) has only one zero ni: 
f{ni) = 0, f{n) 7^ for n 7^ ni. Then, the function F{n) 
has two zeroes: F{0) — and F{ni) — 0. According to 
the Theorem [1] only elements poo, Pnm Pon = Pno can 
have non-zero values in stationary state for the system 
with such NLC. Eq. (|17p implies that stationary values 
of on-zero diagonal elements of the density matrix are 
described by following expressions: 

ni — 1 

POO = E ^'(l"!)' (19) 

fe=0 

where 

g™(|a|) = |are-l"l'/VV^, (20) 

and Pmm 1 - Poo- 

For large enough amplitudes, \a\, of the initial coherent 
state the following estimation is valid: 

Poo < "19^1-1 > 0. (21) 

|a|— foo 

Therefore, the fidelity F — (ni |p| ni) = Pmni of gener- 
ating the Fock state \ni) can be made arbitrarily close 



5 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 




8 10 12 14 

,2 



FIG. 2: Fidelity of generation of the Fock state |ni) by NCL 
starting from the coherent state \a): sohd, dashed, dotted, 
dash-dotted hnes correspond to ni = 1, 3, 5, 7. 



to unity by choosing large enough amphtude |a| of the 
starting coherent state (see Fig. [2]). 



B. Superposition of Fock states |0) and \n) 

As the first example of coherent superpositions that 
can be generated by NCL, we consider states maximally 
close to the state 



|*o„) = 



|0) 



V2 



(22) 



If we require F{n) — F{Q) — and F{m) ^ for 
m ^ 0, n, only elements poo, Pnn, Po« = pJlo (see Eq. ^) 
of the stationary density matrix, p, will be non-zero (The- 
orem [T]). Then, fidelity of generating the state Eq. ((22)) 



F = (*0n I P I *0n) = 



i Re {ponC"'') < i + Ipor, 



(23) 



The maximal possible value of fidelity is F = ^ -I- | pon | ■ 
Positivity of density matrix implies that \pon I < i For 
convenience, we shall characterize coherence of the sta- 
tionary density matrix by the quantity 



COn = 2|pon|, Co„ G [0, 1]. 



(24) 



So, maximization of the coherence co™ leads also to max- 
imality of the fidelity F. 

For the initial coherent state amplitude a = jale*"^/" 
the following equality holds: pon(O) = |pon(0)|e~*'^. Ac- 
cording to Lemma (TJ this holds also for any moment of 
time, t. Therefore, in order to derive conditions for pre- 
serving maximal fidelity by the nonlinear absorption, one 
needs to carry out maximization of the quantity con over 
amplitudes \a\. 



According to Eqs. (d?]), ^ (see also Fig. [rfb)), the 

coherence con equals to 



COn = 2 ^ qm{.\a\)q„i+n{\a\)Tn{T 



(25) 



m— 



where g„i(|a|) are given by Eq. (l20l) . To maximize the 
coherence con one needs to have r„(m) = 1 for m = 
0, n — 1. According to Eq. ^TE\\ . this condition implies 
that 



F{m)=F(m + n) for m = 0, n — 1. 



(26) 



Further maximization of Eq. (j25|) can be carried out 
numerically. Results of numerical calculations are shown 
in Fig. 131 One can see that the best performance of 
the method is achieved for n — 2. Then, according to 
Eq. (23, the coherence is 



co2 = \/2(|a|2-h|a|VV3)e-l"l'. (27) 



The modulus of optimal amplitude of the initial coherent 



state is \aopt\ = y 5 (2 ^ \/3 + V7) ~ 1.2. Elements of 
the final density matrix in the optimal case are poo = 
0.60, P22 = 0.40, IP02I = 0.44 (see Fig. E^b), and the 
coherence is C02 — 0.88. 



C. Superpositions of Fock states |m) and \n) 

The next example is the generation of a state maxi- 
mally close to the state 



V2 



(28) 



Notice that there is rather pronounced difference be- 
tween the current example and the one considered in the 
previous Subsectioir. According to definition of the func- 
tion F{n), F{0) ~ holds for any system. For any am- 
plitude a of the initial coherent state the element poo of 
the initial density matrix has non-zero value. Therefore, 
according to Eq. (jl7p . this density matrix element will 
have non-zero value in the final state. However, its value 
for the ideal state |5'„m) must be equal to zero. 

In order to generate the state sufficiently close to 
|*nm), we require F{m) = F{n) = 0, F{k) ^ for 
k 7^ ri, m, 0. Only the final density matrix elements pij 
with z,j = Q^n,m are non-zero in the stationary state. 
Similarly to the previous example, fidelity of the desired 
state is 



F = I p I 

= i + Re (pn^e'"^) - Poo < ^ 
Again, we define the coherence 



Poo- 



2|p„ 



e [0,1] 



(29) 



(30) 
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FIG. 3: (a) The modulus of the optimal amplitude of the starting coherent state (grey dots) and the maximal achievable 
coherence (black dots) for generation of the state |^'on)- (b) Absolute values of the stationary state density matrix generated 
for the initial coherent state optimal for generation of |^'o2) • 



Further, we use coherence, rather than fidelity, for char- 
acterization of the state preparation quality. 

Phase of the optimal initial coherent state in this case 
is defined by the following equation: 



a e 



I (m— n) 



(31) 



For preserving maximal coherence and for obtaining max- 
imal fidelity of the generated state, we need 



F(k) = F(k + m-n) for fc = 0, 



1 (32) 



(see Eqs. dH]), 

If above conditions are satisfied, the coherence c„„i is 
described by the following expression: 



m — 1 

2 qk{\a\)qk+r. 

k—n 



-n{\a\). 



(33) 



Optimal values of \a\ and maximal achievable co- 
herence Cnm can be found either numerically (dots in 
Fig. |ll[a,b)), or by approximate analytical expressions 
(lines in Fig. |4l^a,b)). 

To obtain an analytical expression for the coherence 
Cnm, we will take into account that (according to central 
limit theorem), for large enough values of jal the follow- 
ing approximation is valid for qfedal) given by Eg. ((20)) : 



1 



■ exp 



(fc-la^^ 



/27r|al ' I 2|al2 
Then, the coherence is approximately equal to 



(34) 



27r a 



^exp 



(fc- |ap)2 + (fc + m-?i- 



(35) 



To maximize the coherence one needs practically the 
same \a\, as required for maximalization of the numera- 
tor in the exponent of Eq. ([55]) . This value is 



'^opt I 



1 

m . 

2 



(36) 



Taking into account Eq. (1361) and approximating sum- 
mation by integration in Eq. (I35L one can derive the 
following expression for the coherence: 



2 erf 



At 



2V2\a, 



exp 



opt I 



8|aopt| 



(37) 



where An = m — n. Fig. |4] shows that obtained approx- 
imate expressions represent results of numerical calcula- 
tions with sufficiently high accuracy. 

The expression Eq. (1571) for the coherence c„m has 
maximal value 0.84 for An^ w 6.4|aoptp. Therefore, the 
method is most suitable for generation of the states with 



(to 



6.4(to 



(38) 



Using the approximation Eq. p4p . one can show that 
the density matrix elements of the optimal stationary 
state are described by the following expressions: 



-'mm — 2 ' 



n = ^erf(2C), 



IPnml =crf(C)e-«', Poo = i{l-crf(2C)} 



(39) 



An 



Examples of density matrices of the states generated 
by the schemes optimized for n = 4, to = 9 and n = 10, 
TO = 17, are shown in Fig. IHc),(d). 
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(a) (b) 




n n 



(c) (d) 




FIG. 4: (a), (b) Optimal amplitude of the starting coherent state (plot (a)) and maximal achievable coherence (plot (b)) for 
generation of the state l^'nm): dots and lines corresponds to numerical optimization and approximate analytical expressions 
(Eq. p6|) . Eq. (I37|) ): black, dark grey, grey, light grey dots and lines correspond to m — n + l,n + 3,n + 5,n + 7. (c), (d) 
Absolute values of elements of stationary state density matrices generated by the schemes optimized for creation of the state 
|*4,9> (plot (c)) and |*io,i7> (plot (d)). 



D. Superpositions of states with equidistant 
photon numbers 

Here we consider examples of functions f{n) (and 
F{n)) with countable sets of zeros. We focus on the case 
of functions with equidistantly distributed zeros: 

f{jN + na) = 0, J =0,1,2,..., (40) 

where N is the distance between neighboring zeros and 
"■Oj < Tio < A^, determines the position of the first zero. 
It should be noted that, according to the definition of 
F{n), one also has F{0) = for any function f{n). 

For such kind of NCL only elements pnm with n,m — 
0, no, no+-/V, nQ+2N, ... remain non-zero in the stationary 
state. 

As in previously discussed examples, for preserving the 
maximal coherence we require that the function F(n) 
should satisfy Eq. (IT51) . This implies 



i.e. function F{n) must be periodic. 

In previous examples zeros of F{n) were giving a clue 
for choosing the amplitude for the initial coherent state. 
It is not so in the current case. Now the amplitude |a| 
will be considered as a free parameter, and the final state 
will be investigated as a function of the amplitude. For 
the sake of simplicity, we assume that the amplitude a is 
real and positive. 

According to Eq. (fT7| , non-zero elements of the station- 
ary density matrix are described by the following expres- 
sion: 

JV-l 

Pnm = ^ qn+k{\a\)qm+k{\a\), (42) 
k=0 

no 

Pon = ^qk{\a\)qn+k{\a\), if no 7^ 0. (43) 

fe=0 



F{n + iV) = F(n), for n = 1, 2, ... (41) 



For large amplitudes of the initial coherent state 
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lap » A'' the approximate expression for g„(|a|), pro- 
vided by Eq. ([M)) . can be used to simplify Eq. (H^ . One 
can show that in this case 



Pn 



Nqn{\a'\)q^{\a'\)-{\ + 0[--^ 



1 



(44) 



where |a'p — |ap — l)/2- Therefore, the stationary 
density matrix can be approximated with 



where 



I* 



("o)| 



(45) 



V. DISCUSSION 

The discussed examples show that systems with ap- 
propriate nonlinear losses can be used for creation of 
Fock states superpositions with sufficiently high fidelity. 
However, for all examples obtained stationary states were 
mixed. It is quite interesting to find the conditions that 
must be satisfied for the final stationary state to be pure 
and to determine whether such conditions can be fulfilled. 

As stated above, Eq. ([T5)) is one of the conditions nec- 
essary for complete coherence preserving. 

Several other conditions are provided by following lem- 



const -5] — 



ijN+no 



(46) 

is the state that can be obtained from a coherent state 
\a') by retaining only the states with photon numbers 
no, no -I- N, uq + 2N, ... On the other hand, the state 
l^^"'') can be considered as a superposition of coherent 
states, distributed on a circle: 



1$ 



(no) 
JV 



) = const 



fc=0 



-2-iTikno / N ^2Tiik / N 



')■ 



For example, for = 2 the states |$2''^) and 1$, 
arc superpositions of states with even and odd numbers 
of photon respectively and correspond to the following 
superpositions of coherent states with opposite phases: 



(47) 



l^f ) ~ W) + !-«') and \<^\'') ~ W) - I - a'). 

To describe "quality" of generation of the superposi- 
tions, it is convenient to consider purity of the final state, 
defined as 



P — Tr p — ^ ^ PnmPmn • 



(48) 



According to Eq. P^ . for large |a|, when presence of 
non-zero element poo can be neglected, this expression 
can be rewritten as 



qn+k^qn+k2 , (49) 



where n = no , no + A^, no -I- 2 A^, .... Approximating sum- 
mation by integration and using Ea. (l34L one can find 
the following expression for the state purity, valid for 
|aP>iV: 




P w I 



1 



24 aP 



(50) 



Therefore, in the limit |q;| — ^ oo the obtained stationary 
state tends to a pure state (Fig. [Sfa)). 

Several examples of density matrix that can be gener- 
ated by the method are shown in Fig. [SJb-d). 



Lemma 3. // all diagonal elements of the initial state 
density matrix are positive, the final stationary state can 
be pure only if the function F{n) has equidistantly dis- 
tributed zeroes. 

Proof. Suppose ni, n'l and n2, n2 are two pairs of succes- 
sive zeroes of the function F(n), and n'j^ — ni < nj — n2. 
Either of n'l, n'2 can be equal to infinity, if the corre- 
sponding zero ni or n2 is the last zero of F{n). 

Then, according to Eq. (IT7|) . the density matrix ele- 
ments have the following values in the stationary state: 



n'^ — ni — 1 

E 

fc=0 



Pni+fe,ni+fe(0): 



E 



(0), 



(51) 



(52) 



Pn2,n2 



n'^—n^ — 1 

Pn2+k,n2+k{^) > 

fe=0 

— ni — 1 

> E 

Pn2+k,n2+k (0), 

fe=0 



(53) 



where we have taken into account positivity of diagonal 
elements of the initial state density matrix. On the other 
hand, the following inequality holds: 

^Pni+k,m+kiO)j (E^"2+fc,n2+fc(0))) > 



> (E \/ Pni+fc.ni+fc(0)/3„2+fc,"2-l-fc(0)) > 



(54) 



> 



E^"i 



+k,n2 + 



fc(0) 



when all the sums are taken over the same range of k. 

Therefore, assumption of nonequal distances n'l — ni 
and n'j — n2 between successive zeroes ni, n'l and n2, n2 
leads to the following inequality: 



Pni,niPn2,n2 |Pni,"2l i 



(55) 



9 




FIG. 5: (a) Purity of stationary states, generated by schemes with the function F{n) having equidistant zeroes: dots (connected 
by solid line) correspond to numerical optimization, dashed lines correspond to the approximate analytical expression (Eq. (|50|) '): 
black dots and lines, grey dots and lines are for A'^ = 2, 3; circles, squares, rhombuses are for no — 0, 1, 2. (b) Matrix elements 
of the stationary state with even photon numbers, generated from coherent state with |Qp — 9 (A^ = 2, no — 0). (c) Matrix 
elements of the stationary state with odd photon numbers, generated from coherent state with \af = 9 (A'^ = 2, no = 1). (d) 
Matrix elements of the stationary state with photon numbers, distributed with interval 3 {\a\'^ — 9, N = 3, no = 1). 



which manifests that the state cannot be pure. Finally, 
the distance between any neighboring zeroes must be 
equal to some constant value, which was denoted by N 
in the last example of Section IIVI □ 



The condition 



Lemma 4. Let ni, n[ and n2, n'2 he two pairs of suc- 
cessive zeroes of the function F{n), n'^ — ni = — n2, 
and let all diagonal elements of the initial state density 
matrix he positive. The final stationary state can he pure 
only if the density matrix of the initial state satisfies the 
following condition: 



Pni+k,ni+k (0) _ (0) 
Pn2+k,n2+k (0) (0) 



for k = 0, n'l — ni — 1. 

(56) 

Proof. Elements Prii,niJ Pni,n2i Pn2.n2 of the stationary 
state density matrix are described by Eqs. ((5T|) . (j52l) and 
the first line of Eq. 



Pni,niPn2,n2 — |Pni,n2l i 



(57) 



which is necessary for the state purity, is satisfied 
only if all parts of Eq. (IM)) are equal to each other. 
Obviously, equality can be achieved only if "vectors" 
(0)} and {pn2+k.n2+k{0)} are coUinear, which 
implies Eq. □ 



On the basis of these lemmas the following theorem 
can be proved. 

Theorem 2. // the function F{n) has at least one zero 
ni, ni > 0, and the initial state is classical, the final 
stationary state will he mixed. 

Proof. Any classical initial state can be represented as a 
mixture of coherent states with positive weights. There- 
fore, if we prove that the stationary state will be mixed 
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for an arbitrary initial coherent state, the stationary state 
will be proven to be mixed for any classical starting state. 

For the initial state being a coherent state |a) the fol- 
lowing statement holds: 

Pn.+Kn.+m ^ |^|2(„.-„,)l!!i±^. (58) 

p„,+fc,„,+fc(0) (ni+/c)! 



Therefore, Eq. (|56)) is not satisfied except for the trivial 
case, when n'^ — ni = 1, F{n) =0. □ 

However, when the amplitude |q;| of the initial coherent 
state is large enough, the generated state can be very 
close to a pure state. Indeed, only the density matrix 
elements Pm^nsCO) '^ith \ni — \a\^\ < |a|, \n2 — < 
\a\ have significantly nonzero values. For ni,n2 3> k 
Eq. ([55)) can be transformed into 



P»i+fc,»i+fc(0) ^ |^|2(n2-ni) f "2 
Pn2+k,n2+k (0) \71i 



a 



2(n2-rii) 



|2(n2-ni) 



1 



112- rii 
ni 



(59) 



1 + 1^ 



Therefore, Eq. (j56p can be satisfied with arbitrarily high 
precision by using high enough amplitudes of the starting 
coherent state. 



VI. CONCLUSIONS 

We have investigated the nonlinear coherent loss as 
a resource for generating non-classical states. We have 
established conditions for generating a prescribed state 
(namely, arbitrarily Fock states and specific superposi- 
tions of them) from the initial coherent state. We have 
highlighted a connection between properties of the Lind- 
blad operator of NCL and the generated state. We have 
demonstrated that the state generated by NCL from the 
initial classical state will always be mixed. However, for 
certain classes of states one can generate almost pure 
states. In some cases it is possible to reach high fidelity 
of the generation by appropriately choosing an ampli- 
tude of the initial coherent state. Fock states belong to 
such a class; one can generate an almost pure Fock state 
for an initial coherent state of a sufficiently high ampli- 
tude. The situation complicates for finite superpositions 
of Fock states. For example, a superposition of two Fock 
states cannot be generated by NCL with the arbitrarily 
high fidelity from the initial coherent state; a strict up- 
per border exists for this case. However, certain infinite 
superpositions can also be generated with an arbitrarily 
high fidelity. For example, one can "comb" the initial co- 
herent state cutting off Fock state components with odd 
number of particles with an arbitrarily high fidelity. 

Finally, it should be noted that though coherent states 
are "nonideal" for generating states via the NCL, they re- 
main optimal initial classical states. Any classical state is 



a mixture of coherent states with positive weights. Quan- 
tum mechanical equations for evolution of the density 
matrix are always linear. Therefore, any final state ob- 
tained from a classical state, is necessarily a mixture of 
final states which are to be obtained from corresponding 
initial coherent states. Thus, no classical initial state can 
lead to purity of the final state, greater than the purity, 
provided by " the best choice" from the possible the initial 
coherent states. 

The authors acknowledge the financial support by the 
BRFFI of Belarus. They are very grateful to V. S. Shch- 
esnovich for helpful discussions. 



Appendix A: Designing the NCL by correlated loss 

To highlight the concept of designing nonlinear loss in 
systems of coupled bosonic modes (which can be realized, 
for example, in Bose-Einstein condensates [Tsl or in op- 
tical fibers [13)1 let us consider a model of A^-l- 1 bosonic 
modes, ai . . . a^+i coupled in the usual linear way to the 
same Markovian reservoir. Thus, the model is described 
by the master equation in the standard Lindblad form: 



dt 



PN+l 



-\Hn+i, Pn+i\+ iJ^N+lpN+l, (Al) 



where Hj^^i is the Hamiltonian describing unitary evo- 
lution of modes described by the annihilation, a^, and 
creation, aj, operators; 7 > is the linear decay rate. 
Here and in the derivation below we define superoper- 
ators Cj on the basis of Lindblad operators Lj in the 
following way: 

CjP — 2LjpL^j — L^jLjp — pL^jLj. 

In Eq. (jAip the Lindblad operator is 

Af-l-l N+1 
Ln+1 = bN+1 = ^ UTV+lj-flj, ^ \UN+Ijf = 1- 

It depends linearly on the annihilation operators aj and, 
therefore, represents a collective mode. One can intro- 
duce a new set of independent bosonic operators, bj, by 
the unitary transformation Uij: 

N+l 

i=i 

where coefficients un+ij are fixed by the interaction be- 
tween the system and the reservoir. The Hamiltonian, 
H]\j^i, can be decomposed in terms of annihilation, In+i, 
and creation, ^jy+ij operators of the collective mode: 



(A2) 



771, n—0 
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where operators Fm,,n depend only on operators bj and 
b] for j < + 1. 

Now we assume that the state of the collective mode 
^AT+i decays to the vacuum very rapidly on the time-scale 
of dynamics prescribed by the Hamiltonian Tf^v+i-Thus, 
an adiabatic elimination of the mode b^+i can be made 
resulting in the following equation for the reduced density 
matrix of modes bj for j = 1 . . . N: 

, . oo . 

J^PN = —^[F,,.,Pm] + J^^^CnPM, (A3) 
n— 1 ^ ' 



where Lindblad operators are L„ = -Fo.n- The mas- 
ter equation (IA3p describes both possible interaction be- 
tween modes bj and their nonlinear losses. Note that 
-Fo.o might include nonlinearities not present in the orig- 
inal Hamiltonian, i?7v+i- 

Even for low order nonlinearities the scheme described 
above is able to lead to the appearance of NCL. In Ref. 
[igj it was shown that the scheme (|Aip for the case of just 
two bosonic modes subject to Kerr nonlincarity leads to 
appearance of NCL with L = an. 
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